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1 Introduction 



For description of low-energy hadron physics phenomenological chiral quark Nambu - Jona- 
Lasinio model [U El E] (see also review [I] and references therein) is successfully applied for 
many years. Mass spectra of scalar, pseudo-scalar, vector and axial vector mesons and their 
low-energy interactions are satisfactorily described in this model [5l[6j[7]. In the simplest case 
of SU(2) x SU (2) chiral symmetry this model contains four arbitrary parameters: current 
mass m of quark doublet u, d (in approximation m u = m^), coupling constant G\ of for- 
quark interaction of the scalar and pseudo-scalar quark currents in the chiral symmetric 
form, coupling constant G2 of four quark interaction of vector and axial-vector currents, 
ultra-violet cut-off parameter A. 

These four parameters allows to describe the pion, the a- meson, vector mesons p and lu, 
axial- vector mesons a±, masses and their main strong decays and weak pion decay constant 
f n . However, only one of these parameters - m - coincides with a parameter of fundamental 
QCD. So the very important problem is to express the rest parameters in terms of QCD 
parameters, e.g. mo and a s (or Aqcd)- But for a long time this problem was not solved in 
a sufficiently satisfactory form. 

In recent works [SJ E] we have succeed in obtaining description of SU (2) x SU (2) NJL 
model using only QCD parameters. This becomes possible due to use of N.N. Bogoliubov 
compensation approach [TUJ [TTJ (application of the approach to QFT problems is described 
in work [12J). As a result a non-local version of NJL model was obtained with uniquely 
defined form-factor. Thus ultra-violet divergences disappear, therefore there is no need of 
introduction of parameter A. Constants G\ and G2 are expressed through m and strong 
constant a s in the non-perturbative region. 

Remind, that application of these results to the sector of scalar and pseudo-scalar mesons 
leads to satisfactory description of 7r and a masses, constant of weak pion decay f n and of 
strong cr — > mi decay. Emphasize, that only parameters m and a s were used. 

It is worth noting, that independent estimate of average non-perturbative value a s was 
obtained in work [T3] (see also [Hj). The same Bogoliubov approach for a study of effective 
non-local three-gluon interaction results in existence of the stable solution for the definite 
form of non-perturbative contributions to running coupling a s (q 2 ). This corresponds to 
average value for the running coupling in the non-perturbative region a s = 0.415 [14J. Taking 
into account this result only one parameter m remains in our disposal. Note, that results 
of works [33], [H] lead to a consistent value of gluon condensate. 

Here we use the version of of non-local NJL model obtained in [5] with the same parame- 
ters mo and a s for calculation of masses and decay widths of vector and axial-vector mesons 
p and a\. Remind that we introduce no new parameters at all. A special attention will be 
paid to value a s = 0.415. 
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2 Compensation equation for effective form-factor 

In the same way as in work [3] we start from the standard Lagrangian of QCD with two 
light quarks and number of colours N = 3 
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In accordance to the approach, application of which to such problems are described in details 
in work [12J, we look for a non-trivial solution of a compensation equation, which is formu- 
lated on the basis of the Bogoliubov procedure add — subtract. Namely let us rewrite the 
initial expression ([1]) in the form 
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2 



\ (7 / - - - - \ 



, } • (jhv^Thv^ + ^757^757^ J + gsihpPAty - (2) 

1 /T 

- j - ^oV F oV) " y • tfA^ - $ *!> $ 1> 
Go f - - - - \ 

G / - - - \ 

- y ■ yPi^ i>ipi> + ^isi^isi^) ■ ( 3 ) 

Here ^ is isotopic doublet, colour summation is performed inside each spinor bi-linear 
combination, F ^ u = d^A u — d u A fl , and e.g. notion G\ ■ ipipipip means non-local vertex in the 
momentum space 

z (2tt) 4 G x F l (pl,p2, p3, pi) 6 {pi + p2 + P 3 + pi) , (4) 

where form-factor i^is introduced, which depends on incoming momenta. The Lagrangian 
contains contribution of both G2 and G3 which are connected correspondingly to isovector 
and isoscalar terms. In the present work we consider compensation equation only for isovector 
four-fermion terms. 

Now we consider the first two lines of the Lagrangian (j2J) as new free Lagrangian Lq, 
and the two last ones as interaction Lagrangian Li nt . Then compensation conditions (see 
again [12], [9] will consist in demand of full connected four-fermion vertices, following from 
Lagrangian Lq, to be zero. This demand gives a set of non-linear equations for form-factors 
Fj. These equations according to terminology of works [TUJ [EE] are called compensation 
equations. 
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In a study of these equations the existence of a perturbative trivial solution (in our 
case Gi = 0) is always evident, but a non-perturbative non-trivial solution may also ex- 
ist. In the present work as well as in previous one we look for an adequate approach, the 
first non-perturbative approximation of which describes the main features of the problem. 
Improvement of a precision of results is to be achieved by corrections to the initial first 
approximation. 

We follow works [8], [9] in definition of the approximation. 

1) In compensation equations we restrict ourselves by terms with loop numbers 0, 1, 2. 

2) In compensation equations we perform a procedure of linearization over form-factor, which 
leads to linear integral equations. It means that in loop terms only one vertex contains the 
form-factor, while other vertices are considered to be point-like. 

3) While evaluating diagrams with point-like vertices diverging integrals appear. Bearing in 
mind that as a result of the study we obtain form-factors decreasing at momentum infinity, 
we use an intermediate regularization by introducing UV cut-off A in the diverging integrals. 
It will be shown that results do not depend on the value of this cut-off. 

4) We use a special approximation for integrals, which is connected with transfer of a quark 
mass from its propagator to the lower limit of momentum integration. Effectively this 
leads to introduction of IR cut-off at the lower limit of integration by Euclidean momentum 
squared q 2 at value m 2 . To justify this prescription let us consider a typical integral to be 
encountered here and perform simple evaluations. Functions which we use here depend on 
variable of the form aq 2 , where a is a parameter having 1/m 2 dimension. 

5) We keep the first two terms of 1/N expansion in equations. 

6) In case of vector vertices there are two Lorentz structures and thus we have generally 
speaking two form-factors instead of one in [H] . However the corresponding set of equations 
has no explicit solution similar to that of [H] and so we proceed in the following way. In our 
approximation we impose simplified kinematic condition that left-side legs of diagrams have 
momenta p and —p, while right-side ones have zero momenta. Now in addition to terms 
proportional to 7 e x 7 e , which we are interested in, terms of the form p x p may be also 
present, Supposing, that the presence of a form-factor connected with the last structure gives 
small corrections we shall transform the initial equation (in diagram form see Fig. 1) to the 
scalar one contracting it with projector of the form 



In the process of the study we have considered also equations obtained with use of projectors 
of more general form, namely 



It becomes clear, that for values d between 1 and 2 the corresponding solutions lead to spread 
of physical values under interest in the range of 5 — 7%, that corresponds accuracy of the 
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method as a whole. So we take the formulated projection procedure as a component of the 
first approximation. 

Now the demand of compensation of full connected four-fermion vertices proportional to 
G 2 multiplied by the vector form-factor leads us to the following equation (see Fig. 1) 
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Here a conversion to Euclidean momentum space is performed, at one-loop level terms 
proportional to iV and 1 are taken into account and for two loops respectively N 2 and 
N. The lower limit of integration is defined by current quark mass m corresponding to 
value «o = 1.925 10~ 8 , which is obtained in the course of consideration of scalar form-factor 
(see [5J [9]). We also use relation 

Gi = ^G 2 ; (8) 

which is derived in the same works. After dividing by G2 that correspond to our intention 
to find a non-trivial solution we integrate by angular variables of four-dimensional space we 
have 
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In view of looking for solutions of Eq. we apply the differential operator 
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to this equation. As a result we obtain a differential equation, which with account of the 
following substitution 

, = e*,P = y a '<l«- 7G '\ do) 

reduces to the following form 
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i.e. it is Meijer equation of the eighth order. Solutions of the equation are represented in 
terms of the Meijer functions [15] with parameters 6j, a*, which we can calculate provided 
Gs is defined. We can naturally admit G$ = G% following rules of NJL model (see jTj). In 
what follows we present confirmation of this assumption. In this case we have 

b x ■= 1.5 ; 6 a = 1 ; b 3 = 0.499991384 ; 6 4 = .500008866 ; 

b 5 = -1.45597130 • 10~ 7 ; b e = ; b 7 = -0.50000003 ; b 8 = -1.0000001 ; 

ai ■= -0.3944464 ; a 2 := 1.9013991 . (12) 

Values of parameters are calculated with account of value mo. 

To obtain a solution of the integral equation we choose four linearly independent solutions 
of Eq. pip decreasing at infinity and form the following linear combination with coefficients 

F(z) = C\G 2 %(^Z 165,64*63,62,61,68,67,66) + ^2 G 28 (% 166,65*63,62,61,68,67,64 ) + 

+ C3 G 2g 1 67,'b4%3 ,62,61,68,66,6s) + ^ 4 ^28 ( z 1 68,J>6%5,64,63,fe2,61,b7) ' (^) 

Coefficients Cj are fixed by boundary conditions, which are obtained in the same way as in 
work pj] 

Q 2 j ^_ g2, \_ Q 2 ^_ Q q\ \ f p f \ 7 _ 

96 2 192 1 384 3 96 2 3 ) tt 4 v^ Xng 
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yF(y)dy = 0; / F (y) = ; / y 3 F (y) = . 
As a result we have 

Ci := 0.3330348455 ; C 2 := 6.254973002 • 10" 8 ; (15) 
C 3 := 3.452159489 • 10" 8 ; C 4 := 2.105889777 • 10" 15 . 

Unlike of scalar case [9] we here do not force the form-factor value at lower integra- 
tion limit to be unity. Using this condition one might try to define ratio of G 2 and G$- 
However assuming equality of these constants we avoid solution of additional complicated 
transcendental equation, but we acquire a criterion of self-consistency of our approach as 
a whole, because calculations show, that changing this ratio in reasonable range we have 
satisfactory results for values of the form-factor at the normalization point. In our case we 
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have F(uq) = 0.96094 and so we consider our assumption to be justified with reasonable 
accuracy. Admissible are values of ratio -§f = X from 1 up to 1.2 as well. For the last 
value F(u ) = 1.098993576. As a matter of fact to fix the ratio one should consider also 
equation for isoscalar vector terms. However this leads to a considerable complication of the 
procedure and so here we only noting, that preliminary estimates show that just for range 
X = 1 — 1-2 values of isoscalar vector form-factor differs from unity not more than by 10%. 
So admitting % = 1 we formulate the ground approximation bearing in mind necessity of 
further corrections. 



3 Wave function of vector states 

We have the non-trivial solution of the compensation equation and thus four-fermion terms 
are excluded from free Lagrangian. There is of course no compensation in interaction 
Lagrangian, which contains these terms with opposite sign. So we can study a problem of 
bound states with account of this four-fermion interaction. Bethe - Salpeter equation for 
vector case in the same approximation as above (see Fig. 2) has the following form . Remind 
that the first approximation corresponds to zero-mass states (in this approximation there is 
the same equation for vector and axial- vector) . 
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Here besides the same kernel as in Eq.Q we take into account also one-gluon exchange 
and one-meson exchange with corresponding constants a s and g 2 /<iit. Note that contribu- 
tions of (pseudo-) scalar mesons here cancel. In Eq. (I16p enters constituent mass m instead 
of current mass in Eq. For parameter m we use results of previous work [5] where it 
was obtained from stability condition for the effective potential. This procedure allows to 
define m corresponding to value of a s . In the same way as in [9] we take values of u — /3 m 4 , 
which correspond to values of a s in the range under study. We perform calculations for 
u = 0.00015, 0.00030, 0.00045. Values of a s , which are slightly corrected in comparison 
with that of work [5] are presented in the summarizing table. 

Differential equation now is the following 

x i z ir^( z wr b ^ Iiz) = - z ( z Tz- ai + 1 )( z i- a ^ 1 h z) - (17) 

where 
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and coefficients hi are roots of the following equation ( Gg — G 2 ) 
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Solution of Eq. (fl7|) decreasing at infinity has the following general form 

^(z) = Cl G*2 8 ^ Ll,'f)2^63,65,f)4,68,67,66^) + ^2^2 8(^161,62^63,64,65,66,67,68 ) + 

+ C*3 C 28 (z lbi'^65,66,63,64,67,68) + ^4 G^ 8 {^Z 1^1^2^3^4 55^7 4,6^8^ + (20) 

riQl( ' |ai,a 2 A 
°5 °2 8 \* 161,62,63,65,66,68,67,64 ^ ■ 

For u = 0.00030 values of parameters hi read 

bx = 1.5; 62 = 1 ; &s = 0.5 ; 64 = 0.763584407 ; (21) 
65 = 0.19323742 ; b 6 = ; 6 7 = -0.7956972342 ; b 8 = -1.161124600 . 

Parameters are the same as before (1121) . Coefficients Cj are defined from the boundary 
conditions 

/•OO /»oo 

^(m 2 ) = l; / *(2/)d3/ = 0; / y^(y)dy = 0; 

J rn 2 J m' 2 

poo poo 

/ y 2 ^(y)dy = 0; y 3 ^(y)dy = 0; (22) 

J m? J m? 

and value g v is given by the iterative procedure being defined by normalization condition in 
one-loop approximation 

Ngl r*(z) 2 F(z) f3 (Gj + 6 G^G 2 ) N 

-az = 1; u = —u; Po = 77 — 7 • (23) 



127T 2 J a Z ' A) 16 7T 4 

Here we introduce into the integral form-factor F (2;) which was obtained in the previous 
section in view to take into account decreasing of interaction for increasing momentum 
variable. 
Ratio 

7 - i < 24) 

Po 754 

gives coefficients for transitions to variable z ~ p A respectfully for vector and scalar sectors. 
Expression for (3 is obtained in works [8], [9] 
With parameters 6, ()2ip we have 

Ci = 1.7465; C 2 = 0.021266; C 3 = 0.00107221; C 4 = 0.00142116; 

C 5 = -0.0000525341 ; g v = 5.00 . (25) 
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4 Results and discussion 

Now we proceed to calculation of observable parameters. Initial estimate of p-meson mass 
is given by expression [7] 

M = ^=; (26) 

where G% is defined from relation (jSJ) and G± is calculated according the method of work [9] 
for chosen value u. Then we introduce one-loop correction to mass squared, which is given 
by expression 

A M 2 = -^7=5 / w w dz ; (27) 
8 tH VP Ju V z 

so that 



M P = ^M 2 + A(M 2 ). (28) 
For values of parameters presented above we have 

M = 974 MeV ; A(M 2 ) = -325908 MeV 2 ; M p = 789MeV . (29) 

We also estimate aimeson mass according to well known relation [7] 

M 2 ai =M p 2 + 6m 2 . (30) 

Coupling constant g p ^2n of p-decay to two n- mesons we find with triangle diagram according 
to the following relation 



3 ~»*.(z?*(JLz) F{z) 



g P ->2K = g s 9v~r^ / dz; (31) 



47T 

where ty s (z) is Bethe Salpeter wave function for scalar states and g s is scalar meson coupling 
according to definition in work [9]. 
The width of p is the following 



^ 2j (M/-4 mi 2 ) 



3/2 



= gp^v H _ " / (32) 

For ai-meson width we consider two channels: a\ — > p7r and ai — > cr7r. The vertex for 
the first decay has the following form (we omit obvious isotopic factor e a t, c ) 

Vp,v{a\ -> P tt) = A) fiV + ^2 ; (33) 

A = 5 / dy, A 2 = / = dy; 

J m 2 y 2 it 2 J m 2 y 2 
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where p, p and q, v are respectfully momentum and Lorentz index for ax-meson and p-meson, 
m is constituent quark mass. 

The second decay is described by the following vertex (isotopic factor <5 a& ) 

V^CLi > G 7r) = g ai ->an (<7 ~ J (34) 



where k and g are respectfully momentum of 7r and a and /i is Lorentz index of a%. Corre- 
sponding partial widths read 

, . (M a 2 1 + M p 2 )(M a 2 1 -M p 2 ) 2 , 2 (M^-igr \ 

4MP5 + A > 16M 2 1 M 2 J ! (35) 

rt 2 r/if 2 - /if 2 i 3 

T-i/ \ yai— «T7r\ ol a J t-i -wis \ | T-i/ \ 

r(ai -> avr) = A87rM 5 ' T <* = F ( ai ~* ^ 7r ) + F ( ai a7r ) ■ 

Here we assume m 2 , << M 2 nrT . 

Using all these expressions we calculate observable quantities for vector and axial- vector 
mesons. Note, that in calculation of widths of decays we substitute calculated masses of 
corresponding mesons. Results are presented in Table. We present there set of calculated 
parameters in dependence on average non-perturbative running coupling a s in range 0.29 - 
0.48. We normalize our calculations by the most precise parameter f n . All other numbers 
in the Table are calculated. We take each column of the Table as a set of the corresponding 
parameters calculated starting from presented there parameters a s and mo- Following the 
general ideology of our approach we consider the last column with a s = 0.415 as the final 
result of our work, bearing in mind, that this value of average a s is obtained in work [T4"] . 
In addition to values of parameters presented in this column we remind parameters of scalar 
sector, which are obtained in the course of performing of previous work [9J. For the same 
a s = 0.415 we have: 

m w = 134MeV; <qq>= -(230 MeVf ; m a = 480 MeV ; T a = 560 MeV . (36) 

Thus the set of parameters seems to be in satisfactory agreement with data. The only 
parameter differing from corresponding experimental value by more than 12% is the mass 
of cti. As a matter of fact a low value for M ai is inherent to other NJL calculations (see, 
e.g. [7]). Presumably in this case there are some additional contributions to be taken into 
account. 

The results being obtained here and in previous work [9] demonstrate that Bogoliubov 
compensation method leads to a reasonable form of effective non-local four-quark interaction 
of the NJL type. As a result this method allows to describe the light meson masses and 
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probabilities of their main decays. Emphasize, that advantages of this approach are absence 
of ultra-violet divergences in quark loop diagrams, inherent in the usual NJL model, and 
the presence of only one arbitrary parameter, namely current quark mass m . Note that in 
spite of this parameter being somewhat larger than its standard value we nevertheless obtain 
reasonable value for constituent quark mass m ~ 260MeV r [9]. Note, that value of mo is 
defined by compensation equation for scalar form-factor in works |8j [9]. Estimates show 
that presence of correction terms in this equation (e.g. the next terms of 1/N expansion) 
changes value of mo significantly, while the observable parameters change only slightly. A 
development of these considerations is the problem for future studies. 

This method in future studies may also be applied to description of electro-weak proper- 
ties of mesons (e.g. pion form-factor, pion polarizability), of tc — it scattering lengths etc. It 
is important problem to expand the approach for chiral U(3) x U(3) symmetry with inclusion 
of the strange quark s. 
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Table 



u 


U.UUUlo 


U.UUUoU 


U.UUU40 


U.UUUoz 


exp /phen 


a s 


n 007 


n /in/i 
U.4U4 


n A SQ 


U.410 


U.410 [14] 




yo 


yo 


yo 


yo 


liipu/c 


9s 


2.66 


2.84 


2.93 


2.86 




m MeV 


21.9 


21.6 


21.2 


21.5 


5-10 


mMeV 


247 


264 


271 


265 


~ 300 


G^ 1/2 MeV 


320 


287 


267 


283 




9v 


4.30 


5.00 


5.52 


5.11 




M p MeV 


713 


785 


830 


791 


771.1±0.9 


Q p 7Y 7T 


4.29 


4.41 


4.36 


4.44 


4.26 


T p MeV 


136 


166 


175 


170 


149.2 ±0.7 


M ai MeV 


935 


1017 


1043 


1018 


1230 ±40 


r oi MeV 


268 


330 


312 


334 


250 - 600 


r(ai -> <rn)/r ai 


0.168 


0.188 


0.201 


0.189 


0.188 ± 0.043 [15] 
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Figure captions 

Fig. 1. Diagram representation of the compensation equation. 

Fig. 2. Diagram representation of Bethe-Salpeter equation for vector bound state. 
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